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Abstract 

We consider solutions to the time-harmonic Maxwell problem in R'^. For such solution we 
provide a rigorous derivation of the asymptotic expansions in the practically interesting situation, 
where a finite number of inhomogeneities of small diameter are imbedded in the entire space. 
Then, we describe the behavior of the electromagnetic energy caused by the presence of these 
inhomogeneities . 

Sur la perturbation de I'energie electromagnetique due a la presence des 
inhomogeneites de petits diametres 

Resume 

Nous considcrons des solutions du probleme harmonique de Maxwell dans R"^. pour une telle 
solution nous obtenons des formules asymptotiques rigoureuses qui sont dues a la presence d'un 
nombre fini d'inhomogeneites avec petit diametre dans I'espace entier. Puis, nous decrivons le 
comportement de I'energie electromagnetique provoque par la presence de ces inhomogeneites. 

Version frangaise abregee 

Dans cette Note, on suppose que dans R'^ on a m inhomogeneites {zj + aBj}"^^, ou a est un 
petit parametre, Bj C R'^ est un ouvert borne et les points {zj}"^^ verifient I'hypothese 
Solent la permeabilite magnetique fia et la permittivite electrique Sa de forme ([2|). 

L'objectif du travail decrit dans cette Note est de voir comment les solutions perturbees 
du probleme ^ se comportent quand un nombre fini d'inhomogeneites {zj + aBj} de petits 
diametres sont introduites dans I'espace entier R'^. Ceci nous amene, egalement, a etudier 
revolution de I'energie electromagnetique correspondante selon cette deformation du milieu de 
propagation. 

Nous commenons notre analyse dans la section 2 en derivant rigoureusement les developpements 
asymptotiques des champs electrique et magnetique ce qui sont uniformement valides dans 
I'espace. Ces formules asymptotiques sont construites par la methode de developpements asymp- 
totiques assorties. Concernant cette methode, le lecteur pent consulter [8l [11] . En domaine 
borne, on peut voir aussi les travaux [2 [31 31 [S]. Le terme d'ordre principale, dans des 
deveoppements asymptotiques analogues (mais pour le cas d'un domaine borne), a ete derive 
par Vogelius et Volkov [IS] et Ammari et al. [5j. Nos formules asymptotiques utilisent des 
tenseurs de polarisation associes a des inhomogeneites electromagnetiques qui semblent etre 
des generalisations normales des tenseurs qui ont ete presentes par Schiffer et Szego jMj et 
completement etudie par beaucoup d'autres auteurs [H [7l [9l [12] . 
Le plan de cette Note est comme suit: 

Dans la section 1, Nous formulons le probleme principale dans ce travail. Dans la section 2, 
supposant que les champs electromagnetiques (Ee,He) sont dans des espaces de Sobolev con- 
venables, nous obtenons a travers la formule de representation integrale de Lippman-Schwinger 
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des developpements asymptotiques uniformes pour les champs electromagnetiques et la densite 
d'energie. Puis, nous formulons une estimation pour I'energie electromagnetique ce qui nous 
apporte aux resultats de stabilite neanmoins la presence de petites inhomogeneites. 

1 Problem formulation 

We suppose that there is a finite number of electromagnetic inclusions in R"^ , each of the form 
Zj + Qf-Bj, where Bj C R'^ is a bounded and C°°-domain containing the origin. The total 
collection of inhomogeneities is: 

-Uj"Li{2j+aSj}. 

The points Zj G R^,j = 1,...,to, which determine the location of the inhomogeneities, are 
assumed to satisfy the following inequality: 

|z,-2,|>co>0,VjV^- (1) 

Let /^o ELud denote the permeability and the permittivity of the free space, > is a given 
frequency and fc — uj^eq^q > is the wave number. We denote by a; = {xi, X2, a^s) the cartesian 
coordinates in R'^ . We shall assume that /ig > and £o > are positive constants. Let > 
and Ej > denote the permeability and the complex permittivity of the j-th inhomogeneity, 
Zj + aBj] these are also assumed to be positive constants. Introduce the piecewise-constant 
magnetic permeability 

Mo(^):=l . ^ (2) 

' \ /ij, X e z-j + aBj, J = 1 ... TO. ^ ' 

The piecewise-constant electric permittivity [x) is defined analogously. If we allow the degen- 
erate case a = 0, then the function iio{x) (resp. £0(2^)) equals the constant /io (resp. £0). 
We suppose now that the collection of inclusions Ba is within an open and fictive subset C R'^ 
and the source current density is at position in R'^ \ f2. 

Let (EQ,Ha) G R'^ X R'^ denote the time-harmonic electromagnetic fields in the presence of 
the electromagnetic inclusions Ba ■ These time-harmonic fields [TOl [13] are the solutions of the 
Maxwell's equations 

V X Eq, = iuj^aiia, in R'^, 

V X Hq = —iujEcMa + Js, in R^, 
X Eq and v x Hq are continuous across d{zj -\- aBj), (3) 

lim \x\ [V X E„ — ik^ x Eq] — and lim \x\ [V x Hq — ik-^ x H„] = 0. 

\x\^oo \x\^oc 

Here v denotes the outward unit normal to d{zj + aBj). We eliminate the magnetic field from 
the above equations by dividing the first equation in ([3|) by fia and taking the curl to obtain 
the following system of equations for Ec,: 



V X ^^^V X Eq - w^£qEq = iwJs, in R^, 

u X Eq is continuous across d{zj + aBj), 

lim |a;| [V x Eq - ik— x Eq] = 0. 



(4) 



Having found the electric field Eq, we then obtain the magnetic field Hq through the formula 

Hq = V X Eq. 

iuj^a 
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The electromagnetic energy is defined by: 

J / ,(£«|Ea(a:,t)P + (^„)-i|H„(x,i)nda;, for t > 0, (5) 

and it is not hard to see that £ e C^(R^). 

The finiteness of the electromagnetic energy requires that both the electric and the magnetic 
field belongs to a space of fields with square integrable curls: 

iJ(curl; R^) := {a e L^(B?f; curl a e L^{B?f}. 

It can be shown that there exists a unique solution E„ e Ci(R+;i^(R^))n(C"(R+;iJ(curl;R^))) 
to the problem ([¥]), and this solution satisfies the following Lippman-Schwinger integral repre- 
sentation formula . 

Lemma 1.1 Let Eq, he the solution of the problem Then, the following integral represen- 
tation formula holds 



m ^ 

E„(x) = Eo(.t) +Y. ( ~ ~ A'")'^' ^ ^(^' • ^o.{x') 

7 = 1 Zj+aBj 



(6) 



+wVo(£i - eo)^(a;, a;') • E„(x')) da;'- 

The 3x3 matrix valued function Q means the Green's function solution to 

V X V X g{x, x') - k^g{x, x') = l3S{x - x'), in R^, 

limi^l^^ |a;| [V x g{x, x') - i/c x g{x, x')] = 0. 

where X3 is the 3x3 identity matrix. In the above notation the curl operator, Vx, acts on 
matrices column by column and V denotes the derivation with respect to the second variable 
x'. 

The proof of Lemma 11.11 follows immediately if we multiply the first equation in (jj]) by g{x, x') -v 
(v S R'^) and if we integrate by parts over the domain which contains the set of inclusions 

2 Asymptotic behavior 

As we said in the introduction our analysis will be dependent on polarization tensors M e R"^^"^. 
We remember that this tensors are defined by 

M{qj/qo;Bj) := / Vv''ix)dx, (7) 

JB-j 

where {qj} (resp. go) denote either the set {sj} or {fij} for 1 < j < m, (resp. denote either Eq 
or /io) and where the functions w', which dependent on the contrast qj /qo, are solutions of the 
following problem, 

V-g(a:)Vt;9(a;) = 0, in R^, 
v'^{x) — X ^ as a; ^ 00. 



go, xeR^\Bj, 
qj, xeBj. 



The function q is given by 

q{x) -- 
The following holds. 

Proposition 2.1 Suppose that and ^ are satisfied and let (Ec,Hc) be the solution of the 
problem (0), then we have 



(i) The electric field satisfy the following uniform expansion 

m 

Ea(a;) = Eo(a;) +a^^[- iuj{fij - no)V' x g{x, Zj) ■ M{p.j/no; Bj)Ho(zj 



(9) 



+cjVo(£j - £0)^^(2;, zj) ■ M{ej/eo; Bj)Eo{zj)] + 0{a^). 
(a) The magnetic field Ha satisfy the following uniform expansion 

m 

H„(x) = Ho(a;) + X! " ^")^' ^ ' ^(ej/eo; 5j)Eo(^j) (10) 

i=i 

-uj^eoifij - fio)g{x, Zj) ■ Mifi-i/fio; Sj)Ho(zj)] + 0{a^). 

To evaluate the influence of the presence of these inhomogeneities with smaU diameters on 
the evolution of the electromagnetic energy defined by ([5]), we shall introduce its density per 
unit of volume which denoted by H^. According to Poynting's theorem, energy density is 
given by: 

-|H„ = V-n„+J, (11) 



dt 



where IIq, is the Poynting vector, 



Mo 



(12) 



The following theorem is concerned with an important result which justifies the behavior of 
the electromagnetic energy when a finite number of inhomogeneities are introduced in the entire 
space. 

Theorem 2.1 Suppose that (QP and ill]) are satisfied and suppose that the inclusion Bj is a 
ball for each j G {1, • • • , m}. Then, the following uniform expansion holds 



-(Ho(x,t)- H,(x,t)) 



m 



£j ~ £0 
£j + 2£o 



. (g;(a;,Zj) •Eo(zj)) 



l/^ioV • [iweoEo X (V' X Q{x,Zj) ■ Eo(zj)) 



k'{g{x,z,)-^o{z,)) X Ho(2 



Mj - Mo 



icj^oJs • (V X Q{x, Zj) ■ Ho(zj)) 



Mi + 2mo 

+ l//ioV • [itJMo(V' X g{x,z.j) ■ Ho(zj)) 
+ fc^Eo X g(a;,z,) •Ho(z,)]l) +0{a'>). 



Proof. 

Recall the following identity: V-(AxB)==B-VxA-A-Vxi?. Relation immediately 
gives 
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V • (n„) = — [H„ • V X E„ - Ea • V X Ha] . 



Mo 



(13) 



Next, under assumption that inclusion Bj is a ball for all j G {1, • • • , m}, it was proved in [7] 
that polarization tensors ([7]) are simplified 



M{pj/no;Bj) 



3mo 
Mj + 2mo 



(14) 



and 
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M(e,/eo; B,) = -^^\B,\I^. (15) 

In other words, the inclusions Bj are symmetric about their centers (balls). Then, according 
to [S| the correction of order four in relation ^ is zero and therefore the remainder is in fact 
0{a^). Using this result and inserting relations ([H]) and (fTS]) into relation the following 
expansions immediately holds: 

E„(x) = Eo(x) +a^J2[~ ^^Mo ^^^i~^°V ,|V' x (?(a;, z,) • Ho(z,) (16) 

^^23(£^Z£2l|s^.|t;(^, • Eo(z,)] + 0(a5), 
In similar fashion we can get an expansion for the magnetic field: 

n^ix) = Ho(x) + a' V [ic^eo^^^^^lS, |V' x ^(a:, z,) ■ Eo(z,) (17) 

-fc'%^|S.|e(-,-.)-Ho(z,)] +0(a^). 
Using relations and p7|) . the following holds 

Tfl 

H„-VxE„ = Ho • V X Eo + a^{^ciHo • V X (V' X ^;(a;,Zj)Ho) 

+ C2H0 • V X {g{x, z,)Eo) + X g(x, z,)Eo) • (V x Eq) 

+ c'^iGix, z,)Ho) • (V X Eo)} + 0{a^), 
where the constants ci, C2, c'^^ and are given by 

^ ^ IT ^ ""^ To — l-^-J'l' 

fc^ iujeq Ej + 2eo 

In the same manner we find the relation for the term Eq, • V x Hq,; therefore relation 
becomes 

m -J 

V-(n„) = V-(no)+a3V — [ciV-((V'xg(a;,z,)Ho)xHo) 

+ C2V • {Gix, 2j)Eo X Ho) + c'lV • (Eq X (V' x g{x, Zj)Eo)) 
+ 4V- (Eo X g{x,z,)Uo)] +0(a5). 
The proof is achieved by expanding the term Jg ■ Eq, at order 5 according to a in (fTTI) . □ 

Based on Theorem 12.11 we can prove the following main result. 

Theorem 2.2 Let T > 0, then there exists some positive constant C such that the following 
estimate as a — > holds uniformly in t ^ [0, T] 

\£^{t) - £o{t)\ < Ca\ 

where the constant C is independent of a and the set of points {zj}JLi provided that assumption 
|IP holds, but this constant C is dependent on \Bj\ and T. 
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